Soft-lubrication effect on the lateral migration of a slightly deformed bubble rising 

near a vertical plane wall 
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Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a 
stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a 
situation with a small clearance c between the bubble interface and the wall. Motivated by the fact 
that experimentally measured migration velocity (Takemura et al. (2002, J. Fluid Mech. 461, 277)) 
is higher than the velocity estimated by the available analytical solution (Magnaudet et al. (2003, 
J. Fluid Mech. 476, 115)) using the Faxen mirror image technique for k(= a/(a + c)) -C 1 (here a is 
the bubble radius), when the clearance parameter e{= c/a) is comparable to or smaller than unit, 
the numerical analysis based on the boundary-fitted finite-difference approach by solving the Stokes 
equation is performed to complement the experiment. To improve the understandings of a role of 
the squeezing flow within the bubble-wall gap, the theoretical analysis based on a soft-lubrication 
approach (Skotheim & Mahadevan (2004, Phys. Rev. Lett. 92, 245509)) is also performed. The 
present analyses demonstrate the migration velocity scales cx Ca e'^Vbi (here, Vbi and Ca denote 
the rising velocity and the capillary number, respectively) in the limit of £ ^ 0. 

PACS numbers: 47.55.D-, 47.15. G-, 47.15.gm, 47.11. Be 



Recent technical progress in generating microbubbles 
[H, including potentials as actuator and sensor, has en- 
hanced the range of appUcations, e.g., additives to reduce 
a turbulence friction drug delivery capsules 0], and 
contrast agents [1]. In many situations, a bubble en- 
counters a boundary wall during its transport process, 
and a hydrodynamic interaction occurs as characterized 
by the inter-scale between the bubble and the wall. In 
practice, it is of primary importance that the bubble un- 
dergoes a repulsive or attractive force in the wall-normal 
direction, which causes a migration and determines 
the bubble distribution, when translating parallel to the 
wall. As the simplest model system, one may raise a 
problem of a spherical bubble rising near a vertical in- 
finite plane wall in a creeping (Stokes) flow. However, 
there is no mechanism to generate the migration force, 
as ensured kinematic reversibility by symmetry of the 
boundary and by linearity in the Stokes equation. In 
fact, the migration force stems from nonlinearities in the 
advective momentum transport 0, 01 and the interfacial 
deformability 0, H] , to break the symmetry. The theo- 
retical success in solving the nontrivial problem of the 
deformation-induced migration of a bubble or drop was 
recently made by Magnaudet et al. 0] using the Faxen 
mirror image technique. However, more recently, an ex- 
periment 0] and a numerical simulation [l3| indicated 
that the theory involves something imperfect since it 
considerably underestimates the migration velocity, de- 
spite of consistencies in terms of the buoyant velocity and 
the deformation. As the most crucial restriction in the 
theory 0], we raise an assumption in the mirror image 
technique that the bubble-wall distance is much longer 



than the bubble radius. In this Letter, to release the 
constraint of the sufficiently wide bubble-wall distance, 
performing experimental and theoretical studies, we ex- 
amine the effect of the distance on the bubble migration. 




FIG. 1: Schematic of a buoyant bubble moving near a verti- 
cal plane wall. The solid outline shows a deformed interface 
with a deflection / from a spherical interface indicated by 
the dotted outline, (a) left panel: coordinates system around 
the bubble. For a bubble radius a and a distance d between 
the bubble centroid and the wall, the geometry is character- 
ized by k(= a/d), which is assumed to be sufficiently smaller 
than unit in the analysis of Ref. 01 . (b) right panel: inter- 
scale coordinates between the bubble surface and the wall and 
scaling relations suitable to a lubrication theory. A clearance 
parameter is defined as e — d/a — 1. 

As schematically illustrated in Fig. [Ija), we deal with 
the migration velocity Vb3{— Vb ■ 63) of a bubble ris- 
ing near a vertical plane wall at a distance d between 
the bubble centroid and the wall. The bubble has an 
equivalent radius a to that of a sphere with the same 
volume. The bubble- wall clearance is given hy c = d — a. 
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The experimental results used here were measured under 
the condition that the Reynolds number 2paVBi/ H (here 
Vsii— • ei), /9, and /i respectively denote the rising 
velocity, the density, and the dynamic viscosity of liquid) 
is unit or less [§]. The pure lateral migration velocities 
induced by the deformation Vbs were calculated from 
them by substituting the velocities induced by the iner- 
tia effects. We also estimated the interfacial deflection 
f{9, (j)) from the sphere by taking the circumference of 
the bubble on the plane X2 = 0. According to the Stokes 
flow theory for the deformation-induced migration [7|, 
the problem is characterized by two parameters, i.e., an 
inverse distance k{= a/d) (or a dimensionless clearance 
e = d/a — 1) and a capillary number Ca = fJ-VRi/j (or a 
Bond number Bo — pa^g/^ as used in Ref. [71] ). Here 7 
and g respectively denote the surface tension and the ac- 
celeration of gravity. Further, as far as Ca <C 1, we may 
use Ca as a perturbation parameter, and reduce the Ca 
dependence in the problem into that Vbi, Vbs/Csl and 
//Ca are dependent only upon k or e. Figure [2] shows the 
migration velocity Vbs away from the wall normalized by 
CaVsi as a function of n. The measured velocity is found 
to be much larger than the analytical solution especially 
for the large k, where the theoretical assumption k ^ 1 
is no longer justified. Such an underestimation of the 
theoretical migration velocity supports the indication in 
the boundary element simulation carefully dealing with a 
drop motion in a shear flow [lo| . The result implies that 
there exists an additional ingredient to generate repulsive 
force for narrow bubble- wall gap, which is not covered by 
the theory Q. 




FIG. 2: Migration velocity Vbs versus k. The solid curve in- 
dicates the analytical solution ,7] Vb3 / {CaVBi) ~ ^40(1+31^/4) 
with the assumption of the sufficiently large distance between 
the bubble centroid and the wall, i.e., k <C 1. The dashed 
curve indicates the experimental results and the circle the 
results obtained by the finite difference simulation discretized 
in the bipolar coordinates system. The inset shows the simu- 
lation result normalized by the analytical solution. 

To clarify the physical mechanism of the repulsive 
force, we numerically address the bubble migration prob- 
lem. We suppose that the bubble quasi-steadily rises 



near an infinite flat plate in an stagnant incompress- 
ible liquid, and both the Reynolds and capillary num- 
bers are sufhciently smaller than unit. Hence, the prob- 
lem is well described by the steady Stokes equation for 
a solenoidal velocity vector and a Laplace law respon- 
sible for the bubble deformation. Free-slip and kinetic 
boundary conditions are imposed on the bubble surface, 
and a no-slip one on the plane wall. The basic equa- 
tions are numerically solved by the second-order finite- 
difference scheme discretized on the bipolar coordinates 
grid, which is boundary-fitted on both the bubble sur- 
face and the wall. Evaluating exact values of the control 
volume and the cell interfacial area, we take care of the 
mass and momentum conservations in a discretized form. 
In a similar manner to Rcf. Q, instead of directly solv- 
ing the flow field with the deformed bubble, we employ 
the Lorenz reciprocal theorem to determine the migra- 
tion force and velocity by coupling two flow fields around 
a spherical bubble translating parallel and perpendicular 
to the wall. To make sure of the numerical stability and 
accuracy, we set the clearance parameter in a range of 
0.1 < £ < 9. 

As shown in the inset of Fig. [21 the ratio of the mi- 
gration velocity becomes close to unit as k approaches 
0, indicating the simulation captures the analytical solu- 
tion for small k [7]. On the other hand, the migration 
velocity for the bubble closer to the wall (i.e., for larger 
k), the discrepancy from the theory is larger. Compared 
with the experiment, the simulation results reveal low 
velocity since it cannot capture unknown factors such as 
small but finite Reynolds number effect, unsteadiness, 
and measurement uncertainty, which may be inevitably 
involved in the experiment. Nevertheless, the simula- 
tion result also indicates the presence of the additional 
narrow-gap repulsive force. 

To demonstrate the narrow gap effect, we investigate 
the near-wall deformation, which causes a change in 
the boundary shape. Introducing an interfacial deflec- 
tion f{9, 4>){= Ca cos cj)), we write the distance 
from the bubble centroid to the interface as a + f. Fig- 
ure [21 shows the angular profile of the scaled defiection 
(Ca)^^2^ Following the analytical solution for n <^ 1 
, it is uniquely arranged as |sin0cos0. In Fig. [SJa), 
the present simulation captures the measured deflection. 
In Fig. [2Ib), the simulation result for the wide bubble- 
wall gap (e — 4) is in agreement with the theory 0]. 
Note that the agreement is confirmed to be better in the 
wider separation. On the other hand, for the narrower 
gap (e = 0.1, 0.5, 0.6), the discrepancy from the the- 
ory is considerably larger especially in the wall neighbor 
{6 ^ 0), that may be related to the larger discrepancy 
in the migration velocity as shown in Fig. [2l Namely, it 
indicates that for small e, the bubble deformation is pref- 
erentially enhanced within the narrow bubble-wall gap, 
and its squeezing effect promotes the bubble migration. 
To shed more light on the role of the inter-scale hydro- 



3 



0.8 
0.4 
0.0 
-0.4 
-0.8 



e=0 6 Simulation 





^ Magnau 


det et al.{2003) 



-;i/4 



1.5 



1.0 



0.5 



0.0 





e 



71/4 



n/1 



Symbols: 
Simulation 

Curves: 

Lubrication Theory 

£=0.1 
--X-- £=0.5 

^ t 



Magnaudet et al.(2003) 



71/4 





3V8 



FIG. 3: Angular profile of the scaled deflection. The 
solid curve indicates the analytical solution [3], — /'"^^'/ft:^ = 
I sin cos 6, with the assumption of k <C 1. (a) upper panel: 
comparison between the experiment (symbols, Ca = 
0.080) and simulation, (b) lower panel: comparison between 
the present simulation, indicated by the symbols, and soft- 
lubrication theory ^ -J^'^^^/k' = ^^^^^^ log(l + 
by the dashed, dotted, and dashed-dotted curves. 



dynamics, a lubrication analysis is performed. It should 
be noted that one can access a relevant physical picture, 
so-called 'soft-lubrication' problem, recently addressed 
by Skotheim & Mahadevan [11], who derived the elastic 
deformation generating the lift force. We prescribe the 
wall-parallel velocity Vbi for ease in completing the anal- 
ysis, and employ standard lubrication assumption, i.e., 
e ^ 1. We also suppose small capillary number Ca ^ 1 
resulting in the small deformation. Further, to make sure 
of the inclination of the near-wall interface to the plane 
wall to be sufficiently small, we adopt an additional con- 
straint e~'^/^Ca <C 1. For the coordinates as illustrated 
in Fig. mb) (see e.g. Ref. 12]), Z = H{R) = 1 + /2 
represents the interface within a portion R 0(1), if 
the deformation is absent. Using e and Ca, we write the 
velocity vector, the pressure, and the deflection / of the 
interface in an expansion form 

^^^^ = ?7W(i?, Z) COS0 + e-^/^Cs£llf\R, Z) + 
Vbi 

(1) 

= Uy{R,Z)sm4> + ..., (2) 



^^^^ = ei/2j7(o)(i?, z) cos(l>+e-^C^UP{R, Z) + . 
Vbi 



«P(!lM = e-3/2p(o) ^ _^ e-3cap(/) (i?) 



/(r, (f)) = e-^''^G&F{R) cos (/), 



(3) 

. . . , 

(4) 
(5) 



whose scaling relations are suitable to all equations to 
appear afterward. Following the spirit of the lubrication 
theory, we solve the governing equations 

1 d{Rur) I du^ ^-l/2dUz _ 

R dR R d(j) dZ ' ^ ' 



dR 



RdS dZ^ 



dp 
dZ 



with the kinetic condition 
i?J7(") - f/i") - R 
-RU, ^g^U, +^j^U^ 2 Qz 

Fdui°^ I f dp f\ ^ ^ 

+ TT— + - + — =0 at Z = H, 

2 dZ 2\dR R 



-0, 
(7) 

(8) 

) 

(9) 



the free-slip boundary conditions 

ac/(°) dM'^ du^J^ Fd^'^ 



dZ 



dZ 



dZ 



2 9Z2 



at Z 



(10) 



and the no-slip boundary condition n = on the plane 
wall Z — 0. Considering that the drag force acting on 
the bubble is deduced O(loge) from the earlier analysis 
for a motion of a rigid sphere |12i] , and the normal stress 
on the bubble surface is dominated by the pressure p ~ 
Q^^-3/2^^ as compared with drUr ~ 0{e~^^'^), Ur/r ^ 
0(£-i/2), u^/r - 0(e~i/2), and d,u, - 0{e-^^^), which 
are related to the viscous stresses, simplifies the Laplace 
law into 



R dR J 



(11) 



We obtain the leading-order pressure and interfacial de- 
flection 



p(0) 



3i? 



F 



31ogF 
5R 



(12) 



To ensure the appearance of the soft-lubrication effect. 
Fig. mb) shows the deflection based on (fT2|) . which is 
in agreement with the simulation results for the narrow 
gap {e = 0.1, 0.5) in the region of the small ^(w e^^^R), 
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indicating the neighborhood of the wall. Therefore, the 
present theory is responsible for the discrepancy of the 
narrow gap experiment and simulation from the available 
theory [7|]. Note that for small e, we confirm that the 
pressure profile for small 9 is also consistent with ([T^. 
From ©J (O and the no-slip condition at Z = 0, we 
write the higher-order quantities as 



UyHR, Z) =A2{R)Z'^ + Ai{R)Z, 



(13) 



p(-^)(i?)= / dR 2^2 (i?), 

J oo 



(15) 



which quantify the deformation-induced migration force. 
The R dependent functions Ai and A2 are determined 
from dH) and Using a solution A2 = -9(4 - 

R^)\ogH/{bOH'^R), we obtain the lateral force Fm to 
cancel the migration velocity and to maintain the wall- 
parallel motion 

Fm lim fiaVei / d(f> U{e^'^R) e^'^R e~^C&P'^^\R) 
Jo Jq 

(16) 



STT/zaCaVBi _2 

it:; £ as e 
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From Bart 13], the drag force acting on the bubble trans- 
lating perpendicular to the plane wall at unit speed is 
Fdc ^ 37r/xa/(2e) as £ ^ 0. Thus, the migration veloc- 
ity asymptotically obeys 

Vb3 = — > — e as e -> 0. (17) 



Fdc 
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To validate the asymptotic solutions, Fig. 2] and the 
inset show the scaled migration velocity and force, re- 
spectively, as a function of the clearance parameter e. 
The simulation results are consistent with two asymp- 
totic behaviors based on the soft-lubrication theory for 
e ^ 1 as well as the mirror image technique for e ^ 1. 
The theories provide the different exponents of the mi- 
gration velocity scaling with respect to £, namely it scales 
Vbs/Vbi oc Ca e^^ for e 1 and Vbs/Vbi oc Ca for 
£ ^ 1. The latter scaling relation is found in this Letter. 

In conclusion, we demonstrated that the former theory 
Q is valid to describe the deformation-induced lift force 
as far as the gap between the bubble and the wall is 
sufficiently wide {k <C 1). On the other hand, for the 
narrow gap case with the clearance parameter e smaller 
than unit, the soft-lubrication effect crucially appears to 
induce the migration velocity asymptotically scalingwise 
Vb3 « Ca 

We thank Shu Takagi for many fruitful discussions. 
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